Two infrared Yang-Mills solutions in stochastic quantization 
and in an effective action formalism 
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Three decades of work on the quantum field equations of pure Yang-Mills theory have distilled 
two families of solutions in Landau gauge. Both coincide for high (Euclidean) momentum with 
known perturbation theory, and both predict an infrared suppressed transverse gluon propagator, 
but whereas the solution known as "scaling" features an infrared power law for the gluon and ghost 
propagators, the "massive" solution rather describes the gluon as a vector boson that features a 
finite Debye screening mass. 

In this work we examine the gauge dependence of these solutions by adopting stochastic quantization. 
What we find, in four dimensions and in a rainbow approximation, is that stochastic quantization 
supports both solutions in Landau gauge but the scaling solution abruptly disappears when the 
parameter controlling the drift force is separated from zero (soft gauge-fixing), recovering only 
the perturbative propagators; the massive solution seems to survive the extension outside Landau 
gauge. These results are consistent with the scaling solution being related to the existence of a 
Gribov horizon, with the massive one being more general. 

We also examine the effective action in Faddeev-Popov quantization that generates the rainbow and 
we find, for a bare vertex approximation, that the the massive-type solutions minimise the quantum 
effective action. 
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I. INTRODUCTION 

The Dyson-Schwinger equations (DSEs) for the gluon 
propagator in Landau gauge have been the object of in- 
tense study in recent years. Two solutions are known in 
several truncation schemes [l| , that we will call "scaling" 
and "massive" . The main difference between them is that 
the first @, Q features an infrared power-law behavior at 
low Euclidean fc 2 while the second produces a propagator 
reminiscent of Yukawa-like theory [3-fToj , with a (gauge- 
dependent) gluon mass. Note however that the gluon's 
Poincare representation is that of a massless vector boson 
with two polarization states, and that any reference here 
to 'mass' should be inter pre ted in the sense of an inverse 
Debye screening length [llj . Note that such confusion 
can be avoided by referring to the massive solution by 
the synonymous term "decoupling" 



12] 



Lattice computations currently seem consistent with 

I 1 T 1 

the second, massive solution [13H18j . but the scaling so- 
lution remains of theoretical interest because of its many 
desirable theoretical properties (the infrared exponents 
can be quite approximately determined for the entire 
tower of Green's functions without truncating the sys- 
tem 0, the Kugo-Ojima confinement scenario [2l[ 
is realized, the non-perturbative BRST quartet [22j ■) 

That these two types of solutions have been found for 
the Yang-Mills system should not have come as a sur- 
prise, as the situation was similar in atomic physics in 
the late 1920's. In the Thomas-Fermi [23| model of the 
atom, a central nucleus and a spherical electron cloud 
act as sources of the electrostatic potential V(r) (akin to 
the ghost propagator in Coulomb gauge), $(x) in appro- 
priately reduced units If the electron energy levels 
are filled according to the Pauli principle in order of in- 
creasing kinetic energy alone, up to the Fermi sphere, one 
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FIG. 1: Solid line (red online): Sommerfeld's power-law so- 
lution of the Thomas-Fermi equation. Dotted line: numeric 
solution to the equation imposing finiteness at the origin. 



obtains the Thomas-Fermi equation 



d 2 <S>(x) 
dx 2 



<S>(xf 



(1) 



If only one boundary condition to this equation is im- 
posed for a neutral atom, namely that the potential van- 
ishes at large distance r a i, one obtains the two well 
known solutions reproduced in Fig. [TJ The solid line in 
the figure is Sommerfeld's solution, $(x) = -±4-, which is 
the unique one with power-law behavior. The dotted line 
represents the usually accepted numeric solution mod- 
elling the atomic charge distribution, which is obtained 
with the additional boundary condition that $(0) = 1 
(other finite boundary conditions lead to additional so- 
lutions that amount to a rescaling of both x and $). 
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Thus, it is natural that modern Yang-Mills Dyson- 
Schwinger equations or Renormalization Group Equa- 
tions, that arc more sophisticated versions of Eq. (QJ, 
also accept solutions that are either finite or power-law 
behaved in the infrared. 

One of the questions we examine in this work is 
whether the existence of the two solutions in Yang-Mills 
theory has any relation to the well-known Gribov ambi- 
guity: the Landau gauge fixing condition d^A^ = does 
not completely fix the gauge, but leaves a discrete num- 
ber of gauge "Gribov" copies J24[- Perhaps, one could 
speculate, the two solutions correspond to different gauge 
representations of the same gauge-invariant information, 
as it is known that the Gribov ambi guity causes a corre- 
sponding ambiguity in the DSE's 2. 



Thus, in the first part of this article we turn to a 
soft gauge-fixing method that avoids the Landau gauge 
condition and the Faddecv-Popov formalism, but pro- 
duces Dyson-Schwinger equations that are quite similar 
in structure to the standard ones. 

The basic observation of stochastic quantization [28| 
is that the weight employed to compute correlators in 
Euclidean field theory 



(A(x)A(y)) = / DAA(x)A{y)e- s ^ YM '> 



[A] 



(2) 



-/3E 



akin to a Boltzmann equilibrium distribution e 
be seen as the end-point e _s ( YM )[ j4 ] = lim r _ j . 00 P(r) of the 
evolution of a stochastic random walk characterized by a 
Fokkcr-Planck equation in a fictitious time parameter r. 
This is 



dP 



d 4 x 



8A a »(x) \5A°(x) 



SP 



K;(x)P 



(3) 



Formally, the force K in this equation is defined to be 
K a Jx) = 

However, Zwanziger showed that the limit r — > oo is 
not well defined in a gauge theory because the stochas- 
tic evolution can run away along a gauge orbit of gauge 
equivalent configurations A(x). To avoid it, one needs 
to add a force term that will pull the evolution towards 
small values of J d 4 x\A(x)\ 2 and thus retain the gauge 
potential at (or near) the closest gauge copy to the origin 
in A-space. This force has to be of dimension 3, to retain 
a hope of renormalizability (the one-loop (3 function has 
indeed been calculated and found standard [29l |30| , and 
the Slavnov- Taylor identities obtained in [3l|), and tan- 
gent to the gauge orbit to avoid disturbing any physical 
observables (expectation values of gauge-invariant oper- 
ators) that may be computed. 

The natural choice is to add a gauge transformation 



(4) 



The parameter a is a real and positive constant that 
controls the relative intensity of the stochastic Yang- 
Mills and the gauge-restoring forces. The gauge is not 



strictly fixed in stochastic quantization (rather, gauge- 
equivalent configurations will be weighted in a smooth 
manner around the origin in gauge space, with much less 
probability for those much farther) except in the limit 
a — > that makes the gauge fixing force dominant and 
fixes the gauge field to Landau gauge. 

The addition of this gauge force, tangent to the gauge 
orbit, makes the method of stochastic quantization geo- 
metrically correct, in the sense that the Gribov problem 
is bypassed. If one then takes the limit r — > oo the path- 
integral weight P, gauge-equivalent to e~ s < YM > , is not 
easy to write down because the new force is not a con- 
servative vector field in A-space, so it can not be written 
as a functional derivative of an action-like functional. At 
least it has recently been shown that the Euclidean solu- 
tion to the time-independent Fokkcr-Planck equation is 
positive and unique (up to normalization) [32|. In con- 
trast, the solution in the Faddccv- Popov theory is highly 
non-unique owing to the nodes of the Faddeev-Popov de- 
terminant at each Gribov horizon (25[. 

However one can undertake the computation of the 
quantum effective action and derive wave equations for 
the Green's functions of the theory. Indeed the Dyson- 
Schwinger equations for the gluon propagator within the 
stochastic formalism have been obtained and described 
in |33| | . 

A peculiarity of this "soft" gauge fixing method is that 
there are no Faddeev-Popov ghosts. Instead one has 
both transverse and longitudinal dressing functions of the 
propagator, with the relative weight a being akin to the 
gauge parameter. 

The gluon propagator is 



S ab D^(k) = I d i x{A a MA b v {x)y i ' 



(5) 



Z T (k 2 ) 



k 2 » 



PUk) + a^P-PUk), (6) 



with 



1 ixv — \ u £2 



P. 



~ k 2 ' 



(7) 



appropriate transverse and longitudinal projectors. Wc 
write down also the "bare" propagator from tree level 
perturbation theory, introducing appropriate renormal- 
ization constants for later use 



D% (fc)-^ 



2 P^v{k) 



(8) 



We will recall the wave equations (in rainbow approxi- 
mation to the DSE of stochastic quantization) for Zl and 
Zt defined in Eq. ((6]) in section |H] The DSE equations 
are then a system of two coupled equations and we per- 
form an infrared analysis in section Hill in order to detect 
possible infrared power-law solutions. This we perform 
in the limit of a — > (Landau gauge) but also for finite 
a. In the first case wc do find a scaling solution, but not 
for finite a. 
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This is corroborated by our numerical solutions to the 
DSE system in section ITVl We find the standard solution 
to the Landau gauge system in Faddcev-Popov quanti- 
zation. Thereafter we solve the corresponding equations 
in Landau gauge in stochastic quantization. Finally we 
proceed beyond Landau gauge and find only the massive 
solution for finite a. 

In section [V] we return to the traditional method of 
Faddcev-Popov quantization, in Landau gauge, where 
both solutions are present. Our goal there is to construct 
an effective action for the traditional rainbow DSE's, and 
evaluate it numerically with the gluon and ghost propaga- 
tor computed in section IIV1 with the idea of identifying 
the absolute minimum among them. We will find that 
the massive solution has least effective action, thus sug- 
gesting a reason why it should be found in lattice gauge 
theory simulations. It appears to us that a constrained 
minimization should be carried out in future work to find 
the scaling solution. What this constraint might be, is at 
present unknown to us, but our work on the stochastic 
quantization method suggests that it may be related to 
the formation of the Gribov horizon. 



Finally, section I VII presents a summary of our results 
and further discussion. 



II. RAINBOW DYSON-SCHWINGER 
EQUATION FOR Z L , Z T 

The Dyson-Schwinger equation for the gluon propaga- 
tor, after projecting via the P T and P L of Eq. (Q, be- 
comes [33[ a system of two coupled equations, see Fig. [5] 



1 _ 7 N 2 

MP) - Zs T Zl9 



u M ( jTTT , otTTL i jTLL 

(2tt)* 



; rr ~Z-iZ n — dZ n Z-\ (J 2 - 

Z L (k 2 ) 3 a ° W 2 



I 111 + 2I IIL + I ILL )(k u k), 

(9) 



a fc l ( T LTT i otLTL , T LLL\ 
(2tt)* 



{I LTT + 2I LTL + I LLL ) (h,k), 



where we employ slightly modified kernels / respect 
to [33j that we now specify; they do depend non-linearly 
on Z^(fc 2 ) and Zr(k 2 ) as usual. To shorten the notation, 
we introduce k 2 = k — k\. Additionally a superscript T or 
L above an index indicates that the index is contracted 
with either of the projectors Pt or Pj, with momentum 
of the corresponding argument; the first of which corre- 
sponds to the external gluon with the second and third 
those gluons internal to the loop. For convenience we 
also write Z3 = Z^Z a and a = Z a a from now on. The six 



kernels read 
I TTT (h,k) 



TTT n M -Z T {k\)Z T {kl) 



h 2 h 2 h 2 

KlZl 2 (fc; -ki, -fe)C(*i. -k) 



(10) 



I TTL (h,k) =■ 



-aZ T (k 2 )Z L (k 2 

h 2 h 2 h 2 



(11) 



KltXik; -ki, -k 2 )G{£ x (h,k2, -k) 



I 1LL {k 11 k) = 



~a 2 Z L (k 2 )Z L (k 2 ) 



1,2 1,2 7,2 

K^X (k; -h, -k 2 )G L x ^ 2X (k u k 2 , -k) 



(12) 



jLTT^^ikDMH) 



h 2 h 2 h 2 
f\j fv2 

(k; -fci, -fel^ft, k 2 , -k) 



(13) 



I LTL (h,k) 



-aZ T (k 2 )Z L (k 2 ) 



h 2 h 2 h 2 



(14) 



KxxX (k; ~ku ~k 2 )GiXx(ki,k 2 , -k) 



I LLL (k u k) =- 



a 2 Z L (k 2 )Z L (k 2 ) 

1,2 1,2 1.1 



(15) 



KtxX -k 2 )G^ 2X (k u k 2l -k) . 

Before application of transverse and longitudinal projec- 
tors, the force that generalises the bare vertex in the 
ordinary DSE's is 



K m3 {k x ;k 2 M = -S(™l 3 (h,k 2 ,k 3 ) (16) 

+ a- 1 A^J 2 ) Al3 (fci;fc 2 ,fc 3 ), 



with 

~^L( fc i' fc 2,fc 3 ) = [(ftiW w]M1 + (cyclic)] , (17) 
4™( fc i5 fc 2,fc 3 ) = [((fc 3 V,W -(2o 3))] . (18) 

These contain the standard three-gluon vertex and the 
Gauge-Fixing term. While the action part is Bose sym- 
metric, the gauge force distinguishes its first index. Fur- 
ther, note that in contrast to Ref. [H[ we have pulled 
out the common factor of (ig) to the front of the DSE 
equation. 

An immediate remark regarding Eq. (| 18[) is that K LTT = 
K TTT = in light of K being proportional to either k 2 
or fc 3 that have zero projection onto the transverse plane. 

Next we specify the dressed vertex G, that essentially 
defines the truncation of the DSE system (we work only 
at one loop in the elementary vacuum polarization in- 
sertion). The G vertex is a sum of a bare three-gluon 
vertex Eq. (fTT|) (that would generate the rainbow ap- 
proximation to the pure Yang-Mills system) and a term 
appropriate to stochastic quantization that incorporates 
the a-dependent drift force, 



GW 2M3 ( fc l, fc 2,fc3) = SjZf£ Ms (kx, fc 2 ,fc 3 ) 

+ r AllAl2Al3 (fci,fc 2 ,fc 3 ) 



(19) 





v\Ow = wvw + 







FIG. 2: Diagrammatical representation of the DSE. The presence of a transverse (longitudinal) projector is represented by a 
spring (wave) in the adjoining gluons. The blobs indicate that the propagator/ vertex is dressed. We have subsumed numerical 
factors and signs into the diagrams. 



where once again, in contrast to Rcf. [331 ] we have pulled 
out the common factor (ig) to the front of the DSE. 

We project each of the indices by trans- 
verse/longitudinal projectors. These quantities are 
anti-symmetric in their three arguments and so 
r ™(^,fc 3 ) = -Tl^ 2 (k 1: k 3 ,k 2 ), and like- 

wise r^j^fci.fca.A*) = -Tl^ikzMM)- These 
stochastic vertices stem from truncated solutions of the 
quantum effective action and are found to be 

t (kl,k 2 ,fo) =0, 

^2 — ^1 



yTTT 



^3 



a{k\ + k\) 



( fc 3) M 



^TLL 



(k 1 ,k2,k 3 ) = — 



k^ &k-^ 



akf 



k\ 



^3 



(20) 
(21) 

(22) 



r 



LLL 



(ki,k 2 ,k 3 ) =— ^ 
a 



x[P T (k 1 )P L (k 3 )]^ 3 -(2^3)y 
(hU (23) 
■ (cyclic) 



k 2 



k 2 



^3 



k-^ ~\~ k>2 
x [P L (k 2 )P L (k 3 )] 



The conventional DSE system of equations equiva- 
lent to Eq. ([9]) in the Faddccv-Popov formalism in Lan- 
dau gauge is long known to have two solutions 0, [34[ . 
The massive solution has a transverse Yukawa-like gluon 
propagator Z(k 2 )/k 2 —> cj(k 2 + m 2 ), and the scaling 
solution instead takes a power-law form Z(k 2 )jk 2 — > 
(A; 2 ) 2k_1 , that, for the (slightly truncation dependent) 
K ~ 0.595 of Faddeev-Popov theory, yields a very sup- 
pressed gluon propagator, while the ghost propagator is 
enhanced G{k 2 )/k 2 -> (A; 2 )"*- 1 . 

A practical way [l2j to numerically select one or the 
other solution is to subtract the second equation of 
Eq. (J9j) at a fixed scale k = to facilitate imposing a 
boundary condition, obtaining 



1 



1 



Z L (k 2 ) Z L (0) 



+ 



(24) 



aZig 



<&k x 



(ij(fci.*Wi(fti,o)) 



(where for brevity we have subsumed the three kernels in 
the integrand into one function Ii ) . Note that subtracting 



at zero is not necessary, rather it increases stability of the 
numerical solution. 

One can now choose Zl(o) ^° ^ e nn ite (in search for 
a massive-like solution) or zero (in search of a scaling 
solution with divergent ghost dressing function) , and try 
to see whether the system does accept both solutions, 
in analogy with the Thomas-Fermi equation, and now in 
the absence of the Gribov ambiguity. 



III. INFRARED ANALYSIS 

To examine the structure of the integral equations near 
a = 0, in search for an extension of the Landau gauge 
power-law scaling solution for the propagators, it is useful 
to classify the various terms according to their power of a. 
Recall that the vertices K and G, in Eqs. p^|) and (|TT))) re- 
spectively, contain a Yang-Mills part together with 
a gauge transformation part (K or T) that can introduce 
a non-trivial a dependence. For brevity, we write the 
kernel product XjJJ 2 (fc, -k 1: -k 2 )Gl^ x (k 1 , k 2 , —k) as 
K\\ A ^AiaTa wnere the ordering dictates the function 
arguments. Then for the transverse gluon we have 



Z T {k 2 )Z T {k 2 ) 

qTTT qTTT 
I \Kl,K) — ; n , n - J/VMWi 



h 2 h 2 h 2 



since K, 



(GF) 
TTT 



I TTL (k u k)=- 



(YM) J (YM) ! 

and Y TTT = 0. For I T ' TL we find 
Z T (k 2 )Z L {k\) ^ 



(25) 



jU2ju2j,2 

r\j fii-j^ luo 



(26) 



?TTL qTLT 



~ ( _ qi r-L nl 

a I ''(YMJ^fYM) 

i ( t^TTL qTLT , j^TTip 
+ ^ A (GF) ,5 (YM) + -"-(GF) 1 - 



qTTL r TLT\ 
"(YM) ) 

ttl^tlt\ 



and finally, 
I TLL (h,k) 



Z L (k 2 )Z L (k 2 

h 2 h 2 h 2 

TfTLL^LLT 
(GF) 



(27) 



aTLL qLLT 
" (YM)°(YM) 



qTLL ^LLT , j^TLL qLLT \ 
*"(YM) + iV (GF)'-'(YM) / | 
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In this last equation we have defined Y LLT = aT LLT in 
order to factor out the explicit factor 5, in Eq. (j2"2")l . 

For the longitudinal gluon propagator, we also include 
the overall factor of a present in Eq. ([9]) to make the 
counting more transparent. We have 



al hll {k x ,k) 



Z T {k 2 )Z T (k 2 ) 

1,21.21,2 

- ( qLTT qTTL _ nLTT f TTl\ 
11 I J (YM)' 3 (YM) J (YM) i J 



(28) 



having used K^q£ — to simplify the expression. 
For I LTL we write 



If a scaling solution should exist, the infrared behavior 
of the dressing functions would respectively be 



Z T (k 2 ) 
Z L {k 2 ) 



■C T (k 2 ) ai 
•C L {k 2 )^ 



(32) 



Substituting into the longitudinal gluon equation (|9|) and 
employing dimensional analysis for small k, with the inte- 
gration dominated by the region of small k\ cx k, we find 
olt + 2(%l = — (4 — d)/2 for (i-spacetimc dimensions, or 
in four dimensions simply the traditional Landau gauge 
result from Faddeev-Popov quantization 

»t = — 2«i . 



aI LTL 



(h,k) = - 



Z T {k 2 )Z L {k 2 ) 

h 2 h 2 h 2 



/ i^ltl^tll\ 



a 2 



(29) 



qLTL qTLL \ 
" D (YM)'-'(YM) J 



~ I tsLTL qTLL qLTL ~y.TLL 
" I A (GF) J (YM) °(YM) i 



once again writing Y TLL = dT TLL . 
Finally, 



For the transverse gluon we must consider the kernels 
of Eqs. ([25 ]) - ([2T]) in Eq. © upon taking the a limit. 
Two of these contributions are given in Ref. [33| 



t t,tt 



(h,k) = 



1 Z T {k\) Z T {k 2 ) 



k\ 



(33) 



Q {YM)TTT n , , s a (YM)TTT n , 

{k,-ki,-k 2 )b{ iX2 \ (*i,«2,-k), 



(34) 



and 



M LLL {k u k) 



Z L (kl)Z L (kl) 
1.21.2 1.2 



a 3 



(30) 



qLLL qLLL \ 

" D (YM) J (YM) J 



~2 I T^LLL qLLL _ qLLL pLLL \ 
a I - fv (GF) (YM) ^(YM) 1 J 



In this last expression we also accounted for factors of a 
in the vertex by writing T LLL = af LLL . 

It is now a simple matter to read off terms that survive 
in the Landau gauge limit, which we investigate in the 
next section. 



A. Landau Gauge 

After substituting the kernels arranged in powers of the 
parameter a (that quantifies the separation from Landau 
gauge) in Eqs. (|2"5j) - (|3"0")) into the Dyson- Schwinger equa- 
tions ([9]) , and reading off the surviving terms in the a — > 
limit, we find for the longitudinal gluon 



1 



1 



g 2 N c 

Z L (k 2 ) Z L (0) ( 2 ^) 4 fc2 



(31) 



k 2 kf-(k-kxY 



k\ (k\ + k 2 ) 



Z T (k 2 )Z L (k 2 ) 



that has already been reported in Ref. |33| where simi- 
larity with the ghost DSE of Faddeev-Popov gauge has 
been noted. 



T,LL 



(ki,k) 



k^e-jh-ky 



Z L (k\)Z L (k 2 ) . (35) 



These two integral kernels are identical to the gluon- and 
ghost-loop kernels in standard Faddccv Popov theory in 
the bare-vertex truncation. However, in Landau gauge 
one contribution that was previously overlooked stems 
from the mixed longitudinal/transverse gluon correction 



t t,tl 



-2 K\ \_ \i{k 



AA 1 A2V t > -h, -k 2 )S XlX2 x(k 1 , fc 2 , -k) (36) 
1 Z T {k\) Z L {k\) 



k\ 



k\ 



This arises from the identification L^ i L A ^ A (fci, ^2, — k)) = 
S? I x T \(ki,k2,—k)) in the limit a — > 0, and hence con- 
tributes additively in G TLT = S TLT + T TLT = 2S TLT 
(this corrects an error in Ref. |3j|, where these terms 
were taken to cancel, which however did not affect the 
calculation reported there). 

Performing the k\ integrals analytically assuming they 
are dominated by the infrared power laws in Eq. (|32l) . 
we find ah = —k — —0.52146. This is consistent with 
an IR enhanced longitudinal gluon propagator Z^/k 2 oc 
(/c 2 ) - 1,52 and an IR suppressed transverse gluon prop- 
agator Zx/k 2 oc (/c 2 ) 04 , similar to the Faddeev-Popov 
ghost-gluon system in Landau gauge but with a slightly 
smaller exponent. 



B. IR analysis for a/0 

We now attempt to find a scaling solution for finite 
To focus only on the important terms, we assume 
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infrared dominance of the gauge- fixing "force" K^ t (A) 
over the Yang-Mills "force" — SSym/^A, in analog to the 
Faddcev-Popov ghost dominance in Ref. [19j . Thus we 
analyse the effect on the integral equation of the gauge- 
fixing force alone, 

K^A) = ar x K gt ^ . (37) 

With this restriction, triply transverse vertices K T,TT 
and r T ' TT vanish, as does K L ' TT . The surviving terms 
in the DSEs are 



«-i J- =1 



L.LL i tL,TL 



1 



T,LL , rT,TL 



(38) 
(39) 



Substituting again the power law ansatze in equation (|32[) 
(with the assumption that no leading-power cancellation 
occurs), we find 



olt = oil , 



(40) 



since independently of which term on the RHS may dom- 
inate we find the same power on the RHS of both equa- 
tions for L and T. By counting powers of k 2 



-a L =d/2-2(l + a L ) 



(41) 



which yields = (4 — d)/6 for Euclidean dimension d. 
Thus 



(XT 


=a L 


= 


for d = 


4 • 


(XT 


=OiL 


= 1/6 


for d = 


3 . 


(XT 


=a L 


= 1/3 


for d = 


2 . 



Thus, in four dimensions we find that the only power- 
law solution is consistent with the perturbative propaga- 
tors proportional to 1/k 2 , the dressing not carrying an 
anomalous infrared dimension [50|. 

Other studies concerning scaling behaviour of the 
Yang-Mills system in gauges different to Landau gauge 
have been performed, such as ghost-antighost symmetric 
gauges HU , maximally Abelian gauges [H, H3] and inter- 
polating gauges [38| . Similarly, linear covariant gauges 
have also been investigated on the lattice [3^, H(| . Typi- 
cally these introduce an additional scale which dominates 
in the IR. However, note that the naive scaling analysis 
does not preclude the existence of a mixed scaling regime, 
in which the scaling behaviour is only manifest at inter- 
mediate, rather than infrared, momenta. An example of 
this scenario is QED in 2+1 dimensions 
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Upon numerically solving the complete system of equa- 
tions, see Eq. ©, matching the infrared and the ultra- 
violet behavior with a computer, we expect to find both 
the massive and scaling solutions for a = 0, but perhaps 
only the massive solution for a ^ 0. We thus turn to the 
computer in the next section. 



IV. NUMERICAL SOLUTION OF THE DSE 
EQUATIONS IN STOCHASTIC QUANTIZATION 

As a warm-up we first show the solution to the conven- 
tional system of the coupled one-loop gluon-ghost DSE 
equations that has been discussed at length elsewhere 
recently [T2|. Note here that in contrast to recent stud- 
ies, we employ the bare vertex approximation since we 
will later investigate the effective action for this Faddeev- 
Popov system. In contrast to Ref. [42| which also em- 
ployed the bare vertices approximation, we will not em- 
ploy an angular approximation. 

For the sake of clarity, we briefly overview the solution 
method. For the scaling solution, one assumes the ex- 
istence of power-law behaviour for the ghost and gluon 
dressing functions: 



G(p 2 ) = Ap- K Z{p 2 )=Bp 2 



(42) 



with the precise value of k dependent upon the assump- 
tion of ghost dominance and the Taylor condition Z\ = 1 . 
Through the running coupling associated with the ghost- 
gluon vertex one surmises the existence of an infrared 
fixed point in the coupling. This gives a stringent rela- 
tionship between the coefficients A and B in Eq. (|42[) . 
Typically, one chooses a value for g 2 and A, with B then 
determined through knowledge of the IR fixed point . The 
renormalization constants for the ghost and gluon prop- 
agators are removed by subtraction, in favour of renor- 
malization conditions for Z(/j, 2 ) and G(/i 2 ). 

The ghost equation must be subtracted at /x 2 = zero 
to choose the boundary condition, as in Eq. (|2~4"|) . To 
search for the scaling solution we impose the condition 
1/G(0) = in accordance with Kugo-Ojima. The final 
condition for Z(s 2 ) is then surmised by smooth matching 
of the numerical solutions to the IR. Ultimately, fi is then 
arbitrary and one scales the solutions in the momentum 
variable to match say the physical value of the running 
coupling at some scale. 

The resulting ghost and gluon dressing function for 
both massive and scaling solutions are represented in 
Fig. [31 

In the case of Stochastic QCD we follow the procedure 
as closely as possible. The longitudinal gluon is sub- 
tracted at zero momenta, whilst the transverse gluon is 
subtracted at some large perturbative scale. In the Lan- 
dau gauge limit one similarly has access to an IR fixed 
point. By employing bare- vertices we lose Multiplicative 
Renormalizability and so, it is less trivial to relate the 
choice of the renormalization condition Z(s) to the mo- 
mentum scale. 

To address the massive solutions, one merely replaces 
the boundary condition for G(0) = finite. Note that we 
could perform the subtraction of the ghost and gluon 
equations at the same large perturbative momentum, 
obtaining both massive and scaling solutions. In this 
case fine-tuning of the renormalization condition dictates 
which type of solution is selected. 
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p [arbitrary units] p [arbitrary units] 



FIG. 3: Numerical solution to the system of coupled Faddeev-Popov ghost-gluon equations at one loop (rainbow approximation 
with bare vertices, in Landau gauge). Both scaling and massive solutions are represented. Left graph: ghost dressing function. 
Right graph: gluon propagator dressing function. 



One technicality of this particular system of equations 
is the appearance of quadratic divergences in the trans- 
verse equation for the gluon propagator. These are elim- 
inated by merely employing an additional subtraction in 
the infrared. 

In Fig. 0] we then plot the analogous quantities in the 
framework of stochastic quantization in Landau gauge, 
the longitudinal and transverse dressing functions, set- 
ting a — > (Landau gauge). The numerical result for 
the scaling solution fulfills our expectations based on 
the analytical infrared scaling study, with the correct 
k value that we obtained analytically. The massive so- 
lution seems a simple modification of that obtained in 
Faddeev-Popov quantization, which comes as no surprise 
as the two systems of equations in Landau gauge, whether 
stochastic or Faddeev-Popov, are very similar. 

Finally we proceed to finite (but very small) a in Fig. [5] 

Our numeric results confirm that the scaling solution 
with a suppressed transverse gluon is a feature of stochas- 
tic QCD only in Landau gauge, and that a massive-type 
solution can however be found even with soft gauge fix- 
ing. This is in agreement with extant lattice data [43ll44| ]. 



V. EFFECTIVE ACTION FOR THE 
FADDEEV-POPOV SYSTEM 

In this section we return to Faddeev-Popov quantiza- 
tion of Yang-Mills theory and consider it well established 
that functional methods find two solutions to the wave 
equations. We plot these in a simple truncation in Fig. [3] 

It is thus puzzling that lattice gauge theory reproduces 
only one of them. A natural avenue of investigation is 
to consider an effective quantum action from which the 
functional (be it DSE or Renormalization Group equa- 
tions) are derived, and study its value for the different 
solutions. Perhaps lattice gauge theory is picking up the 



absolute minimum of the effective action and the second, 
scaling solution, can only be found with a constrained 
minimization. Recently, a study of the renormalization 
group [45J singled out the 'massive' solutions as those sta- 
ble in the infrared. Here, we wish to investigate whether 
any of these solutions is a local minimum of the effective 
action, or rather if they are all saddle points. 

A convenient starting point is the two-particle irre- 
ducible (2PI) quantum effective action of Yang-Mills the- 
ory to one- loop. Propagators in the effective action are 
the ones in the fully interacting vacuum, but the ver- 
tices are taken from perturbation theory. Thus, taking a 
functional derivative respect to the explicit propagators 
produces the rainbow Dyson- Schwinger equations to one 
loop without need to worry about implicit dependences 
of the vertices. It is for this reason that we employed 
the bare-vertex approximation in the study above. The 
effective action reads [46| 

T (D, G) = *-Tr log D- 1 + % -Tr (D^D) 

- iTrlogG- 1 - iTt (G^G) + T 2 [D,G] , 

(43) 

where the non-trace part involving higher than 2-point 
functions is given as 

T 2 [A G] = +^DXl ~ y 2 DG 2 V ( f )2 , (44) 

and represented diagrammatically in Fig. [5] The vertices 
here are V03, that is the three-gluon vertex in perturba- 
tion theory (at this order identical with the Lagrangian- 
level vertex in Eq. (JTFJ S™f 3WJ (ki,k 2 , k 3 )) and the bare 

ghost-gluon vertex Vq£ . 

The full gluon and ghost propagators in Landau gauge 
are symbolically represented by D and G, and their per- 
turbativc counterparts by Dq and Go- 

If one takes a functional derivative of Eq. (|4"3")l respect 
to the propagators 5T/5D = 0, ST/SG = 0, the rainbow 
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FIG. 5: Numerical solution to the system of coupled longitudinal-transverse gluon equations at one loop (rainbow approximation 
with bare vertices, in stochastic quantization). The scaling solution now has disappeared (only the perturbative propagator 
comes out of the analysis), while the massive solution still exists. Left graph: longitudinal dressing function Zl- Right graph: 
transverse gluon dressing function Zt- 




•oooooooc* 



FIG. 6: Non-trace part (involving dressed three-point func- 
tion) of the effective action that generates the rainbow Dyson- 
Schwinger equations of Yang-Mills theory by seeking the sta- 
tionary point under variation of G, Z (diagrammatically, cut- 
ting through each line). 



one-loop Dyson-Schwingcr equations for D and G follow. 
For example, let us derive the abstract expression for the 



ghost equation 
ST(D,G) 



SG 



This gives 



iG n 



G 







G-^G^+g^DGV^ 2 



l -g 2 D-2GV t )2 



(45) 
(46) 



and likewise, for the gluon equation 



D- 1 = Dn 1 + ^-D 2 V r 2 



(47) 



With the propagators parametrized in terms of the 
conventional gluon and ghost dressing functions 



G(k) = 



k 2 



G(fc 2 ) 
fc 2 



(48) 
(49) 



the effective action is a functional of Z(k 2 ) and G(k 2 ) 
that has zero variation at the solution of the Dyson- 
Schwingcr equations. The reduction of the effective ac- 
tion to a simple form that can be handled by a computer 
given Z, G is shown in appendix |A"1 The outcome is, up 
to constant contributions, 



T =8V 



d 4 k 



log (G~ 1 (k 2 ))+G(k 2 ) 



(50) 



2 9 



3 
~ 2 
d A q 



(27T) 
d A q 



[log(Z-\k 2 ))+Z(k 2 )] 
4 K 1 (k,q)G(k 2 )G(q 2 )Z((k-q) 2 ) 



(2tt)< 



K 2 (k,q)Z(k 2 )Z(q 2 )Z((k-q) 2 ) 



with kernels 



Ki(k,q) 
K 2 (k,q) 



k 2 q 2 — k ■ q 2 
k 2 q 2 r A 

k ■ q 2 — k 2 q 2 
k 4 q A r 4 



[3r 2 (k 2 



(51) 

2k 2 q 2 + k-q 2 }, 



where r = k — q. 

A very appealing way of visualizing the effective action 
is to take a curve in function space that passes by the two 
solutions of the DSE equations in Landau gauge (massive 
and scaling). The OY-axis of Fig.[7]is the effective action 
Ta(D, G) for cutoff A, whereas the OX axis represents an 
arbitrary interpolation parameter a that varies between 
and 1. 

So that we can show also the propagators on the per- 
turbative vacuum state Go, Z$, we choose a parabola 
through all three functions (we actually employ an inter- 
polation of their logarithm since the functions are very 
dissimilar, and all three are positive) 

log Z[a] = (2a - 1) [(a - 1) log Z Q + a log Z 2 ] 

-4a(a-l)logZi, (52) 

log G[a] = (2a - 1) [(a - 1) log G„ + a log G 2 ] 

-4a(a-l)logGi, (53) 

obtaining, presumably a local maximum at Zq and then 
an absolute minimum at either of (Z\, G\) or (Z 2 , G 2 ). 

We thus proceed to express the effective action V in 
terms of generic ghost and gluon dressing functions, Z 
and G, with the eye on evaluating it over the parabola 
passing by the three relevant choices. We will drop all 
constant contributions to the action, as we are interested 
in the relative ordering of the solutions and not the value 
of the action itself, which is of little relevance. 
Our numeric results are shown in Fig. [JJ 
We separately show the contribution from the kinetic 
terms (Ff roc , denoted with a dashed line) that features 
a maximum of the action (along this slice) at the per- 
turbativc propagator (on the left of the plot at a = 0). 




-5e-05 



FIG. 7: The solid line depicts the 2PI cut-off quantum effec- 
tive action (units GeV 4 per unit four-volume) that generates 
rainbow-DSE equations, plotted on a slice through the func- 
tion space of ghost and gluon Euclidean momentum propa- 
gators. For a = the propagator is bare, for a = 1/2 its 
dressing functions yield the massive solution, and for a = 1 
we arrive to the scaling solution. We choose a cut-off that 
emphasises the difference in the infrared. 



This is natural as this function satisfies the free Dyson- 
Schwingcr equation, thus it must be an extremum of the 
free effective action. 

That is no more the case for the full computation in- 
cluding interactions (solid line) where the perturbative 
propagator plays no special role. However a clear min- 
imum can be seen very near the massive solution at 
a = 1/2, that seems to be the extremum of the inter- 
acting action. The scaling solution is not an extremum 
of the unconstrained minimization. Note that one should 
not attribute meaning to the precise location of the min- 
imum in Fig. [7J as this depends upon the interpolating 
function used. It is merely indicative of which class of 
solutions minimises the action. 

Finally, the figure separately depicts F 2 , the interact- 
ing part of the effective action at two loops, as defined in 
Eq. (|31]), that is the difference between the total and the 
free actions. In conclusion, from an analysis of the effec- 
tive action in bare- vertex approximation, it appears that 
the massive solution to the Dyson-Schwinger equations is 
naturally found in lattice gauge theory due to its lesser 
action. Whether the scaling solution can be found by an 
adequately constrained minimization, as well as a more 
extensive exploration of the effective action, deserves an- 
other publication. 



VI. SUMMARY 

We have addressed the contemporary topic of the two 
solutions in the Dyson-Schwinger equations of pure Yang- 
Mills theory, after pointing out that finding both finite 
and infrared-scaling solution has really been with us since 
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the time of Thomas' analysis of the atom, since Sommer- 
feld actually found a power-law solution to the Thomas- 
Fermi equation for the electrostatic potential there in ad- 
dition to the well-known finite solution. 

In Yang-Mills theory which of these two solutions is 
found depends on the chosen gauge. We have employed 
Landau gauge to show within a standard rainbow trun- 
cation scheme, in conventional Faddeev-Popov quanti- 
zation, that the massive solution has less effective ac- 
tion than the scaling solution, and thus we conjecture 
that this is the reason why the massive (or decoupling) 
solution to the gluon DSE is so easily found in lattice 
gauge theory. Whether this remains to be the case for 
the dressed system is yet to be seen, but we expect the 
qualitative picture to remain the same. 

For our main contribution we have adopted the view- 
point of stochastic quantization, where a gauge fixing is 
not forced, but instead the system adopts a statistical 
distribution where, in each gauge trajectory, it tries to 
relax to the gauge copy closest to the origin. The control 
parameter of this thermodynamic gauge force, a, when 
set to 0, allows to guarantee the Landau gauge condition 
in stochastic quantization. Within this Landau gauge, we 
find results similar to the Faddeev-Popov method, with 
both scaling and massive solutions. 

If this condition is lifted, so that Landau gauge is not 
fixed, by allowing a to be finite, the scaling solution to 
the Dyson-Schwingcr equation abruptly disappears. The 
massive family of solutions remains as the only one for 
finite a. Thus we have another piece of evidence that 
might suggest that the scaling solution is related to the 
Gribov horizon forming in the curvilinear Landau gauge. 

In conclusion, we believe, from this and other works, 
that there are indeed two classes of solutions to the com- 
plete wave equations of Yang-Mills theory, as found in 
various truncations of the Dyson-Schwinger and the Ex- 
act Renormalization-Group equations. We also conjec- 
ture that lattice gauge theory is finding the solution with 
least effective action in an unconstrained minimization, 
which is a massive-like solution, and that the scaling so- 
lution probably needs to be obtained with a minimiza- 
tion that takes into account additional constraints. Our 
study of the two solutions with the stochastic quantiza- 
tion method finds that the scaling solution is a property 
of Landau-gauge fixing, and disappears if configurations 
not in Landau gauge arc allowed to contribute to the path 
integral. This gauge dependence is in line with other find- 
ings in the literature |48|. We conjecture here that the 
scaling solution is related to the formation of the Gribov 
horizon and that perhaps this is the direction to search 
for an appropriate boundary condition that restricts the 
lattice configurations in order to also find a possible scal- 
ing solution in Monte Carlo simulations of Yang-Mills 
theory. 
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Appendix A: Evaluation of the Effective Action for 
Faddeev-Popov quantization in Landau gauge 

1. Evaluation of the free action 

We first deal in this subsection with the terms that are 
of zero order in the strong coupling g. 

Consider the two terms of Eq. (jl3"]l that include the 
logarithm of an inverse propagator D av or G. Note that 
log A is best interpreted in the diagonal basis 



log 



A x 



A, i 



logAi 



log A n 



(Al) 



so that, for a constant times the identity, log (c<5 afc ) = 
(logc) • 5 ab . Likewise, log (f(a)S ab ) = (log(/(a))) S ab . 

We can easily calculate the traced logarithm of the bare 
ghost propagator Gq 1 = Gq lab (x, y) = i\3 x 5(x - y)8 ab 
to be 

iTrlogGo- 1 = I d A x8 aa {\ogiU x ) S(x - y) 

a •* x=V 

(A2) 

The constant contribution log(ix) = logi + log a; = log a; 
can be neglected, and we represent the delta function by 



8(x - y) = 



d A k 
(2^j 



I ex P [—ik(x — y)] , 



(A3) 



that takes an additional factor i in Euclidean space be- 
cause d 4 k -> id 4 k E . Since £ Q 6 aa = 8, / d 4 x = V with 
V the space-time volume we find 

r (fife 

- tTrlog Gq 1 = 8-V j log(fc 2 ) + constant (A4) 

J (27r) 

Now we dress the ghost propagator, k 2 — > qjjjt) (de- 
noted in co-ordinate space by a tilde, D x — > D x ) and 
obtain 

-^'logG-^/l^log^. (A5) 

Next we address the traced logarithm of the gluon 
propagator, ^TrlogZ?" 1 . Starting again with its bare 
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counterpart and ignoring constant contributions, 
^Trlogi^ 1 = ^Trlog - W x S ab S(x - y)] 



4 r 

4iJ2S^ d 4 xlog [(S^a - d^) x ] S(x - y) 

K=l J 



y=x 



(A6) 



Once again we replace the delta function by Eq. (|A3|) and 
set x = y to leave 



-TrloglV 1 



(A7) 



We write log [-<5^/c 2 + k^k v ] = log(-l) + log(/c 2 (5^(fc)) 
in terms of the projector transverse to k, 5^ u (k) = 
Sfj,v — k^ky = log(k 2 )8 fJlI/ (k) (it is not difficult to convince 
oneself that this relation is valid in the _Lfc subspace writ- 
ing the expression in coordinates. The constant log(— 1) 
is dropped). Thus 

4 

^TrlogZ?- 1 = ^-vJ d 4 kY,$»Ak)logk 2 



-24 



■ V 



(2tt)' 



log k 



(A8) 



Proceeding to the dressed gluon similar to that of the 
ghost yields, 



^TrlogD" 1 = -12V / -^rlog 
2 6 J (2tt) 4 * 



A: 2 



Z(k 2 



(A9) 



Next we address the quadratic terms in the free action, 
-zTr (G^G)) = 

-»£ / / d i xd 4 yGo lab (x,y)G ba (y,x) = 



ah 



Si 



d xd y^- 



U x ff d 4 k d 4 q 



(2tt) 4 (2^) 4 



-i(k-q)(x-y) 



(A10) 



Note the 'dressed' d'Alembertian to represent the inverse 
propagator in position space. Because of translational 
invariance we make a change of variables with unit Ja- 
cobian, z — x — y, z' = ^^ L 1 with the integral over z' 
yielding the integration volume V. Thus after simplifica- 
tion with the Euclidean conventions we have 



iTr (Gq 1 G) = 8V 



d 4 k 



G(fc 2 ). 



(All) 



Likewise we construct 



ab fiv 



§EE // d 4 xd 4 y(5 ab {5, v U-d,d v )J{x-y) 



x S ba (d^n-d^ ^ 1 S(x - y)j . (A12) 
Passing to Fourier momentum space and simplifying 

: (1 + 3- Z(k 2 )) . (A13) 



— Tr [Dq D) = —V , 

2 K ' 2 J ( 27r ) 4 



Collecting the free parts of the effective action together 
we find 



/^JlogfjL + G(k 2 ) (A14) 



(2tt) 4 1 °\G(k 2 ) 



3, / k 2 \ l + 3Z(k 2 ) 
2 l0 H^)J 2 



in which we can discard the constant contribution to the 
effective action that results from log(fc 2 ). 

It is easy to check that this free effective action, upon 
minimization, demands the G and Z ghost and gluon 
dressing functions to take their tree-level value 



SG{q 2 ) 
SZ(f) 



const. (-G- 1 + 1) 



Gfree{q 2 ) — 1 , 

/3 3 
™t.^-Z q -- 

Zfree(<l 2 ) = 1 • 



(A15) 



(A16) 



2. g part of the effective action 

We now turn to the interacting part of the action, 
Fig. |6] In rainbow approximation it is given by the sum 
of two terms, -^g 2 D 3 V^ 3 and - \g 2 DG 2 V^ h)2 '. Let us 
start by computing the ghost loop 



2 nr f2 T/ (gfc)2 



2^ 03 



2 9 



d 4 k 



d 4 q 



iG(k 2 )iG(q 2 )iZ({k-qY 



k 2 



q 2 (k - qf 



8^ 



(A17) 

(if abc k^) (if acb q u ) 
{k-q)*(k-qy \ 
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where we employ the conventions of Rcf . [46j . Summing 
Euclidean and color indices, 



12g J 



(A18) 



d 4 k d 4 q k 2 q 2 ~{k-qy 
(27r) 4 (27r) 4 k 2 q 2 {k-q) 4 

G(k 2 )G(q 2 )z({k-qf) . 



Because of invariance under rotations, the external k- 
integral can be reduced to 27r 3 / °° dkk 3 , and the re- 
maining invariance under two rotations of the inter- 
nal g-integral (whose integrand depends on cos 9 = 
k ■ q since k has already been fixed) simplifies to 
47T J °° dqq 3 d9 sin 2 8 which leaves a total of three inte- 
grations to eventually be performed on a computer. 
Whereas the momentum integrals display a quadratic 



divergence that need to be regulated consistently with 
the DSE, the angular integral is convergent at end points 
due to the sin 9 factor. 

We reduce the gluon loop following the same proce- 
dure. 

After simplifying (with the aid of FORM Hz|), 
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9 2 D 3 V \ = ^ I 



4 ,^Z(k 2 )Z(q 2 )Z(r 2 ) 



k 4 q 4 r 4 



(A19) 



(k ■ q 2 - k 2 q 2 ) (3r 2 (fc 2 + q 2 ) + 2k 2 q 2 + (k ■ qf) X , 



where r — k — q Again, this is an eight-dimensional in- 
tegral featuring quadratic divergences in each k, q; it is 
reducible to 3D as before. 



A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D 65 [18 
(2002) 025012 >ep-ph/0107078| . 

L. von Smekal^ ~K~. Hauck and R. Alkofer, Annals [19 
Phys. 267, 1 (1 998) [Erratum-ibid. 269, 182 (1998)] 
[hep-ph/9707327] . 

L. von Smekal, R. Alkofer and A. Hauck, Phys. Rev. Lett. [20 
79 (1997) 3591 [hep-ph/9705"242] . 

J. E. Mandula and M. Ogilvie, Phys. Lett. B 185, 127 [21 
(1987). 

J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). [22 
D. Dudal, S. P. Sorella, N. Vandersickel and H. Ver- 
schelde, Phys. Rev. D 77 (2008) 071501 [arXiv: 071 1.44961 [23 
[hep-th]]. 

D. Dudal, J. A. Gracey, S. P. Sorella, N. Vandersickel [24 
and H. Verscheld e, Phys. Rev. D 78, 065047 (2008) [25 
[arXiv:0806.4348l [hep-th]] . 

A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. [26 

D 78 (2008) 025010 |arXiv:0802. 18701 [hep-ph]] . 

P. Boucaud, J. P. Leroy, A. Le Yaouanc, J. Micheli, [27 

0. Pene and J. Rodri guez-Quintero, JHEP 0806 (2008) [28 
099 [arXiv:0803.216Tl [hep-ph]] . [29 
A. C. Aguilar, D. Binosi and J. Papavassiliou, JHEP 
1201, 050 (2012) [arXiv:1108.5989l [hep-ph]]. [30 
R. Alkofer, PoS FACES Q CD (2010) 030 [31 
[arXiv:1102.3166l [hep-th]] . 

C. S. Fischer, A. Maas and J. M. Pawlowski, Annals [32 
Phys. 324, 2408 (2009) [arXiv:0810. 19871 [hep-ph]]. 
A. Maas, A Cucchieri and T Mendes, Braz. J. Phys. 37, [33 
219 (2007) [hep-lat/0610006l . 

A. Cuc chieri and T. Me ndes, Phys. Rev. D 81, 016005 [34 

(2010) [arXiv:0904.4033l [hep-lat]]. 

E. -M. Ilgenfritz, C. Menz, M. Muller-Preussker, 

A. Schi ller and A. Stern beck, Phys. Rev. D 83, 054506 [35 

(2011) [arXiv:1010.5"T20l [hep-lat]]. 

1. L. Bogolubsky, E. -M. Ilgenfritz, M. Muller- 
Preussker and A. Sternbeck, PoS LAT 2009, 237 (2009) [36 
[arXiv:0912. 22491 [hep-lat]] . 

O. Oli veira and P. Bicu do, J. Phys. G G 38, 045003 [37 
(2011) [arXiv:1002.4l5T1 [hep-lat]]. 



O. Oliveira, P. J. Silva and P. Bicudo, PoS FACESQCD 

, 009 (2010) |arXiv:1101. 59831 [hep-lat]]. 

R. Alkofer, C. S. Fischer and F. J. Llanes-Estrada, Phys. 

Lett. B 611 (200 5) 279 [Erratum-ibid. 670 (2009) 460] 

[hep-th/0412330] . 

C. S. Fischer and J. M. Pawlo wski, Phys. Rev. D 75 
(2007) 025012 [ hep-th/0609009] . 

P. Watson and R. Alkofer, Phys. Rev. Lett. 86 (2001) 



5239 [hep-ph/0102332 
N. Alkofer and R . Alkofer, Phys. Lett. B 702, 158 (2011) 
[arXiv:1102. 27531 [hep-th]] . 

Llewelyn H. Thomas, Proc. Camb. Phil. Soc 23, 542 
(1927). 

D. Zwanziger, Nucl. Phys. B 321, 591 (1989). 

D. Zwanziger, Phys. Rev. D 65, 094039 (2002) 

[hep-th/0109"224] . 

A. Maas, Phys. Lett. B 689, 107 (2010) |arXiv:0907.5T85l 
[hep-lat]]. 

A. Maas, larXiv: 1106. 39421 [hep-ph] . 

G. Parisi and Y. -s. Wu, Sci. Sin. 24, 483 (1981). 

A. Munoz Sudupe and R. F. Alvarez-Estrada, Phys. Lett. 

B 166, 186 (1986). 

A. Munoz Sudupe, Phys. Lett. B 167, 221 (1986). 

R. F. Alvarez-Estrada and A. Munoz Sudupe, Phys. Rev. 

D 37, 2340 (1988). 

N. Vandersickel and D. Zwanziger . larXiv: 1 202 . 149 1 II [hep- 
th]. 

D. Zwanziger, Phys. Rev. D 67, 105001 (2003) 
[hep-th/02 06053 . 

P. Boucaud, J-P. Leroy, A. L. Yaouanc, J. Micheli, 
O. Pene and J. Rodri guez-Quintero, JHEP 0806 (2008) 
012 |arXiv:0801.272"Tl [hep-ph]]. 

R. Alkofer, C. S. Fischer, H. Reinhardt and 
L. von Smekal Phys. Rev. D 68 (2003) 045003 
[hep-th/0304134] . 

M. Q. Huber, K. Schwenzer and R. Alkofer, Eur. Phys. 
J. C 68 (2010) 581 |arXiv:0904. 18731 [hep-th]]. 
R. Alkofer, M. Q. Huber, V. Mader and A. Windisch, PoS 
QCD -TNT-II (2011) 003 [arXiv: 11 12.61731 [hep-th]]. 



13 



[38] C. S. Fi scher and D. Zwanziger, Phys. Rev. D 72 (2005) 

054005 [hep-ph/0504244] . 
[39] A. Cucchieri, A. Maas and T. Mendes, Comput. Phys. 

Commun. 180 (2009) 215 arXiv:0806.3l24l [hep-lat]]. 
[40] A. Cucchieri, T. Mendes, G. M. Nakamura and 

E. M. S. Sant os, PoS FACESQCD , 026 (2010) 

[arXiv:1102. 52331 [hep-lat]] . 
[41] C. S. Fischer, R. Alkofer, T. Dahm and P. Maris, Phys. 

Rev. D 70 (2004) 073007 [hep-ph/0407T04] . 
[42] D. Atki nson and J. C. R. Bloch, Phys. Rev. D 58 (1998) 

094036 |hep-ph/9712459] . 
[43] J. M. Pawlowski, D. Spielmann and I. -O. Stamatescu, 

Nucl. Phys. B 830, 291 (2010) |arXiv:0911.492T1 [hep- 
lat]]. 

[44] H. Aiso, J. Fromm, M. Fukuda, T. Iwamiya, A. Naka- 
mura, M. Stingl and M. Yoshida, Nucl. Phys. Proc. 
Suppl. 53, 570 (1997). 

[45] A. Weber, larXiv: 1205.04911 [hep-th] . 

[46] J. Berges, Phys. Rev. D 70, 105010 (2004) 



[hep-ph/0401172l . 

[47] J. Kuipers, T. Ueda, J. A. M. Vermaseren and J. Vollinga, 
larXiv: 1203.65431 [cs.SC] . 

[48] A. Maas, Ph ys. Rev. D 79, 014505 (2009) 
[arXiv:0808.3047l [hep-lat]] . 

[49] For completeness we note that x = r/a, with a = 
ao(97r 2 /(128Z)) 1/ '' i , in terms of the nuclear charge Z 
and Bohr's radius ao, and $ = rV(r)/(Ze) rescales 
the Coulomb potential due to the spherically symmet- 
ric (though not pointlike) charge distribution. 

[50] As a side remark we comment that the disappearance of 
the scaling solution in four dimensions in the setting of 
stochastic quantization is perhaps expected because of 
the analysis of 0] since the scaling solution seems to be 
a feature of ghost-antighost symmetric Faddeev-Popov 
quantization, and such symmetry is not implemented in 
the stochastic scheme. 



